Event-plane flow analysis without non-flow efl"ects 
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The event-plane method, which is widely used to analyze anisotropic flow in nucleus- nucleus col- 
lisions, is known to be biased by nonflow effects, especially at high pt- Various methods (cumulants, 
Lee- Yang zeroes) have been proposed to eliminate nonflow effects, but their implementation is te- 
dious, which has limited their application so far. In this paper, we show that the Lee- Yang-zeroes 
method can be recast in a form similar to the standard event-plane analysis. Nonffow correlations 
are eliminated by using the information from the length of the flow vector, in addition to the event- 
plane angle. This opens the way to improved analyses of elliptic flow and azimuthally-sensitive 
observables at RHIC and LHC. 

PACS numbers: 25.75.Ld,25.75.Gz,05.70.Fh 



I. INTRODUCTION 

Since elliptic flow has been seen at RHIC it has 

become a crucial observable for our understanding of the 
system created in heavy-ion collisions. Anisotropic flow 
is most often analyzed using the event- plane method Q ■ 
These analyses are plagued by systematic errors due to 
nonflow effects Q. Nonflow effects are particularly large 
at high Pt [3], where they are likely to originate from 
jet-like (hard) correlations; they are expected to be even 
larger at the LHC. The purpose of this paper is to show 
that nonflow effects can be avoided at the expense of a 
slight modification of the event-plane method. 

Anisotropic flow of an outgoing particle of a given type, 
in a given phase-space window, is defined as its azimuthal 
correlation with the reaction plane Q 



Vn = (cos(n((/) - <^r))) 



(1) 



where n is an integer {vi is directed flow, V2 is elliptic 
flow) , (j) denotes the azimuth of the particle under study, 
<i>i{ the azimuth of the reaction plane, and angular brack- 
ets denote an average over particles and events. Since '^r 
is not known experimentally, w„ cannot be measured di- 
rectly. 

The most commonly used method to estimate Vn is the 
event-plane method [J|. In each event, one constructs 
an estimate of the reaction plane <I>/j, the "event plane" 
Hi- The anisotropic flow coefficients are then esti- 
mated as 



D„{EP} 



1 

'r 



;cos(n(0-*fl))) 



(2) 



where R = (cos(n(5'fl — $/?))) is the event-plane reso- 
lution, which corrects for the difference between r and 
i^R. This resolution is determined in each central! ty class 
through a standard procedure [9']. 

The analogy between Eq. ([2]) with Eq. ([T]) makes the 
method rather intuitive, but its practical implementation 
has a few subtleties: 



• One must remove autocorrelations: the particle un- 
der study should not be used in defining the event 
plane, otherwise there is a trivial correlation be- 
tween and ^ fi Q. This means in practice that 
one must keep track of which particles have been 
used in defining the event plane, so as to remove 
them if necessary. 

• Nonflow correlations between the particle under 
study and the event plane must be avoided. This 
cannot be done in a systematic way, but rapidity 
gaps are believed to largely suppress nonflow ef- 
fects [3,[i3|. 

• Event-plane flattening procedures must be imple- 
mented to correct for azimuthal asymmetries of the 
detector acceptance 

A systematic way of eliminating nonflow effects is to 
use improved methods such as cumulants pd| or Lee- 
Yang zeroes [r2||. Cumulants have been used at SPS [l3| 
and RHIC p, [l4ll. Lee- Yang zeroes have been imple- 
mented at SIS^ and at RHIC fTi|. The reason why 
these methods have not been more widely used is proba- 
bly that they are thought to be tedious. 

In this paper, we show that the method of flow anal- 
ysis based on Lee- Yang zeroes can be recast in a form 
similar to the event-plane method. Specifically, Eq. ^ 
is replaced with 



?;„{LYZ} = {WRCOB{n{cp - 



(3) 



where ^"7? is the same as in Eq. ([2]), and Wr is a new 
quantity, the event weight, which depends on the length 
of the flow vector. The advantage of this improved event- 
plane method over the standard event-plane method is 
that both autocorrelations and nonflow effects are auto- 
matically removed. 

The paper is organized as follows. In Sec. [TTl we de- 
scribe the method for a detector with perfect azimuthal 
symmetry, and we explain why it automatically gets rid 
of autocorrelations and nonflow correlations, in contrast 
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to the standard event-plane method. Readers interested 
in applying the method should jump to Appendix [XJ 
which describes the recommended practical implemen- 
tation, taking into account anisotropics in the detector 
acceptance. In Sec. lIIIl we present results of Monte-Carlo 
simulations. Lee- Yang zeroes are compared to 2- and 4- 
particle cumulants. Sec. IIVI concludes with a discussion 
of where the method should be applicable, and of its lim- 
itations. 



II. DESCRIPTION OF THE METHOD 
A. The flow vector 

The first step of the flow analysis is to evaluate, for 
each event, the flow vector of the event. It is a two- 
dimensional vector Q = {Qx,Qy) defined as 

M 



Qx — Qcos{n'^ii) = ^^Wj cos{n(f)j) 

M 

Qy = Qs\TL{n'^ = "^^Wj sin{n4>j), (4) 
i=i 

where the sum runs over all detected particles, M is the 
observed multiplicity of the event, 4>j are the azimuthal 
angles of the particles measured with respect to a fixed 
direction in the laboratory. The coefficients Wj in Eq. @ 
are weights depending on transverse momentum, particle 
mass and rapidity. The best weight, which minimizes the 
statistical error (or, equivalently, maximizes the resolu- 
tion) is Vn itself, Wj{pT, y) cx Vn{pT, y) [I3|- ^ reasonable 
choice for elliptic flow at RHIC (and probably LHC) is 
w = pt- 

Lee- Yang zeroes use the projection of the flow vector 
onto a fixed, arbitrary direction making an angle nO with 
respect to the a;-axis. We denote this projection by Q^: 

= Qx cos n9 + QySmnO = Qcos{n{^ii -6)). (5) 



B. Integrated flow 

The first step of the analysis is to measure the inte- 
grated flow Vn , defined as the average value of the pro- 
jection of Q onto the true reaction plane: 



K = (Qcos(n(«'fl-$fl))) 



(6) 



where angular brackets denote an average over events. 
We define the complex-valued function: 



1 



(7) 



The modulus \G^{r) \ can be plotted as a function of r for 
positive r. In the case where collective flow is present, 



this modulus has a sharp minimum, compatible with 
within statistical errors \\^. This is the Lee- Yang zero, 
which is determined by finding the first minimum numer- 
ically. We denote its value by . The estimate of the 
integrated flow is then defined by 



Joi 



(8) 



where joi — 2.40483 is the first zero of the spherical 
Bessel function Jo (a;). 

Please note that the above procedure only makes use 
of the projection of the flow vector onto an arbitrary di- 
rection 9. For a perfect detector, azimuthal symmetry 
ensures that is independent of 0, up to statistical er- 
rors. In practice, however, it is recommended to repeat 
the analysis for several values of 9 (see Appendix 



C. Event weight 



Once is determined, the event weight is defined by 



WR = -J^{r'Q), 



(9) 



where Ji{x) is the spherical Bessel function of first order, 
Q is the length of the flow vector, and C is a normal- 
ization constant which is the same for all events in the 
centrality class. This constant is determined by requiring 
that the projection of the flow vector Q onto the reaction 
plane corresponds to the integrated flow Vn, defined by 
Eq. ([6]). This gives the condition 



{WrQ) = 
Inserting Eqs. ([8|) and ([9]), we obtain 



C=—{r'QJ,{r'Q)) 
Joi 



(10) 



(11) 



where angular brackets denote an average over events. 

This normalization constant can be shown to be related 
to the resolution parameter x used in the standard event- 
plane analysis (see Appendix I A 31 for a precise definition): 



Ji(joi 



(12) 



Fig. [T] displays the weight Wr as a function of the 
length of the flow vector Q, together with the proba- 
bility distribution of Q, for two values of the resolution 
parameter. The normalized distribution of Q is Q 



dN 
dQ 



■ exp 



For X ^ 1, this distribution is a narrow peak centered at 
Q^Vn- The weight defined by Eqs. (HD and ^ is then 
close to 1 for all events. If x is smaller, the distribution of 
Q is correspondingly wider, and Wr is negative for some 
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FIG. 1: (color online) Shaded area: probability distribution 
of Q, Eq. with Vn = 0.0625 (see Sec.|lITll. Open circles: 

histograms of the distribution of Q obtained in the Monte- 
Carlo simulation of Sec. IIIII following the procedure detailed 
in Appendix [X] Solid curve: weight Wr defined by Eqs. (|9]) 
and H12p . Stars: weights obtained in Sec. IIIII Top: x = 1-5, 
corresponding to a resolution R = 0.86 in the standard anal- 
ysis (see Eq. (HI). Bottom: x = !> corresponding to a reso- 
lution R = 0.71. This is the typical value for a semi-central 
Au-Au collision at RHIC analyzed by the STAR TPC 



events. These negative weights are required in order to 
subtract nonflow effects, as will be explained in Sec. Ill Dl 
On the other hand, they also subtract part of the flow. In 
order to compensate for this effect, the global normaliza- 
tion of the weight increases when x decreases (compare 
the curves in the top and bottom panels in Fig.[T]). This 
qualitatively explains the x dependence in Eq. (H^). 

The weight ([S]) vanishes linearly at Q = 0. This is 
physically intuitive. The idea of the flow vector is that 
by summing over all particles, one increases the relative 
weight of collective flow over individual, random motion 
of the particles. If the flow vector is small in an event, 
it means that the random motion hides the collective 
motion in this particular event, which is therefore of little 
use for the flow analysis. 



D. Autocorrelations and nonflow effects 

We now explain why the method is insensitive to au- 
tocorrelations and nonflow effects. We consider the situ- 
ation where there is collective flow in the system, so that 



Vn is non-zero for most of the particles. Next, we assume 
that the particle under study does not take part in collec- 
tive flow, i.e., it has w„ — 0; on the other hand, it may be 
correlated with a few other particles (for instance if they 
belong to the same jet). Such nonflow correlations usu- 
ally result in u„{EP} ^ (generally > for intra-jet cor- 
relations). By contrast, our improved estimate i;„{LYZ} 
vanishes, as we now show. 

We separate the flow vector, Eq. ([4]), into a flow part 
Qf , involving particles which contribute to collective 
flow, and a non-flow part Qnf, involving the particle 
under interest (autocorrelations) and the few particles 
which are correlated to it, but do not contribute to col- 
lective flow: 

Q = Qf + Qnf. (14) 

We further assume that the flow and the nonflow part 
are uncorrelated. 

Eqs. ^ and ^ deflne our estimate of Vn as 

v„{LYZ} = 1 ( Ji(/Q) cos(n(0 - ^-fl))) . (15) 

We rewrite the Bessel function as an integral over an- 
gles [H 

Ji(r''Q)cos(n(</>-^'^)) = / |^e- « cos(n(,/>-(?)), 

Jo 

(16) 

with defined by Eq. (O. Since the flow vector appears 
in an exponential, the flow and nonflow contributions to 
Eq. p3|) can be written as a product of two independent 
factors: 

.„{LYZ} = -A £^ I (e-^^.^^ cos(n(0 - 9)) . 

(17) 

We then use again the assumption that the flow and non- 
flow part are uncorrelated to write 

(e-^0^) = (e-^«?)(e-^Q.^.). (18) 

Now, (e''^"^'^^ = by definition of the Lee- Yang zero r^, 

up to statistical fiuctuations. This means that one of the 
two factors in the right-hand side of Eq. ([T8|) vanishes. 
Since it is the flow which produces the zero, it means that 

e""'''5t ^ = 0. Inserting into Eq. we flnd 

t;„{LYZ} = 0, (19) 

up to statistical fluctuations. This completes our proof 
that the method is not biased by autocorrelations and 
nonflow effects. 



III. SIMULATIONS 

= 28 000 events were simulated with a Monte-Carlo 
program dubbed GeVSim|20|]. In GeVSim the V2 and 
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the particle yield as function of transverse momentum 
and pseudorapidity can be set with a user-defined pa- 
rameterization. For these simulations events were gen- 
erated using a linear dependence of V2{pt) in the range 
0-2 GeV/c, above 2 GcV/c the V2{pt) was set constant. 
The average elliptic flow is {V2) = 0.0625. We then re- 
constructed V2(j>t) from the simulated events using sev- 
eral methods: the Lee- Yang-zeroes method described in 
Appendix |X1 2- and 4-particle cumulants The corre- 
sponding estimates of W2 are denoted by W2{LYZ}, W2{2} 
and W2{4}, respectively. W2{2} is generally close to V2 
from the traditional event-plane method; both arc biased 
by nonflow effects. On the other hand, W2{4} is expected 
to be close to i;2{LYZ}, with the bias from nonflow ef- 
fects subtracted. The weight Wj in Eq. ([4]) was chosen 
identically 1/M for all particles, with M the event multi- 
plicity, so that the integrated flow Vn defined by Eq. ^ 
coincides with the average elliptic flow, i.e., Vn = 0.0625. 
The analysis was repeated twice by varying the multi- 
plicity M used in the flow analysis: the values 256 and 
576 were used, so as to achieve a resolution of x =1 a-nd 

1.5. m 

Fig.[2]shows the generated (input) V2{pt) together with 
the reconstructed V2{pt) using cumulants and Lee- Yang 
zeroes for x — ^- The upper panel shows the results in 
the case where all correlations are due to flow. In this 
case, all three methods yield the correct V2{pt) and {V2) 
within statistical uncertainties (see Tabic which are 
twice larger for W2{4} and W2{LYZ} than for W2{2} (see 

Sec. Ell- 
in the lower panel, simulations are shown which in- 
clude nonflow effects. Nonflow correlations arc intro- 
duced by using each input track twice, roughly imitat- 
ing the effect of resonance decays or track splitting in 
a detector. Experiments at RHIC have shown Q that 
nonflow effects are larger at high-pt (probably due to jet- 
like correlations), and a realistic simulation of nonflow 
effects should take into account this pt dependence. Our 
simplified implementation, which does not, is not realis- 
tic. It is merely an illustration of the impact of nonflow 
effects on the flow analysis. Fig. [2] shows that due to 
nonflow effects, the method based on two-particle cumu- 
lants (?;2{2}) overestimates the average elliptic flow {v2)- 
The error on the average elliptic flow is larger than 20% 
(see Table m right column). The transverse-momentum 
dependence of W2(pt) is also not correct, with an excess 
at low Pt by 0.03. By contrast, the results from 4-particle 
cumulants (w2{4}) and Lee- Yang zeroes (w2{LYZ}) are, 
within statistical uncertainties, in agreement with the 
true generated flow distribution. This shows that the 
method presented in this paper is able to get rid of non- 
flow effects. 



IV. DISCUSSION 

Two effects limit the accuracy of flow analyses at 
high energy: nonflow effects and eccentricity fluctua- 




FIG. 2: Differential elliptic flow U2(pt) reconstructed using 
different methods: the line is the input V2. The upper, lower 
panel shows the obtained V2{pt) for the different methods 
without and with nonflow present, respectively. 



TABLE I: Value of the average elliptic flow {V2) reconstructed 
using different methods with and without nonflow effects in 
the simulated data. The input value is {v2) = 0.0625. 



Method 


Flow only 


Flow-f nonflow 


«2{2} 


0.0626 ± 0.0003 


0.0764 ± 0.0004 


W2{4} 


0.0624 ± 0.0005 


0.0627 ±0.0007 


i;2{LYZ} 


0.0626 ± 0.0005 


0.0629 ± 0.0007 



tions |2l|, 123] ■ The method presented in this paper 
is an improved event-plane method, which eliminates 
the first source of uncertainty, nonflow effects. It has 
been recently argued [13] that cumulants (and there- 
fore Lee- Yang zeroes, which corresponds to the limit of 
large-order cumulants) also eliminate eccentricity fluctu- 
ations However, a recent detailed study [l^l 
shows that even with cumulants, there may remain large 
effects of fluctuations in central collisions and/or small 
systems. This issue deserves more detailed investigations. 

Letting aside the question of fluctuations, we now dis- 
cuss which method of flow analysis should be used, de- 
pending on the situation. There are three main classes 
of methods: the standard event-plane method [3], four- 
particle cumulants [ll| , and the Lee- Yang-zeroes method 
presented in this paper. When the standard event-plane 
method is used, nonflow effects and eccentricity fluctua- 
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tions are generally the main sources of uncertainty on ti„, 
and they dominate over statistical errors. The magnitude 
of this uncertainty is at least 10% at RHIC in semi-central 
collisions; it is larger for more central or more peripheral 
collisions, and also larger at highpt. Unless statistical er- 
rors are of comparable magnitude as errors from nonflow 
effects, cumulants or Lee- Yang zeroes should be preferred 
over the standard method. 

The advantages of Lee- Yang zeroes over 4-particle cu- 
mulants are: 1) They are easier to implement. 2) They 
further reduce the error from nonflow effects. 3) Effects 
of azimuthal asymmetries in the detector acceptance are 
much smaller. 4) The statistical error is slightly smaller 
if the resolution parameter x > 1. For x = 0-8, the er- 
ror is only 35% larger with Lee- Yang zeroes than with 
4-particle cumulants (and 4 times larger than with the 
event-plane method). 

Our recommendation is that Lee- Yang zeroes should 
be used as soon as x > 0.8. For small values of x, typi- 
cally X < 0-6, statistical errors on Lee- Yang zeroes blow 
up exponentially, which rules out the method; the sta- 
tistical error on 4-particle cumulants also increases but 
more mildly, and their validity extends down to lower 
values of the resolution if very large event statistics is 
available. 

A limitation of the present method is that it does not 
apply to mixed harmonics: this means that it cannot 
be used to measure vi and V4 at RHIC and LHC using 
the event plane from elliptic flow [25|. Note thatih can 
in principle be measured using Lee- Yang zeroes [26| us- 
ing the "product" generating function, but this method 
cannot be recast in the form of an improved event-plane 
method. Higher harmonics such as W4 also have a sen- 
sitivity to autocorrelations and nonflow effects, which 
is significantly reduced by using the product generating 
function \u\ . 

Although we have only explained how to analyze the 
anisotropic flow of individual particles, it is straightfor- 
ward to extend the method to azimuthally dependent 
correlations [13, HI] . The only complication is that the 
azimuthal distribution of particle pairs generally involves 
sine terms (29j . in addition to the cosine terms of Eq. ([T]). 
These terms are simply obtained by replacing cos with 
sin in Eq. ©. 

In conclusion, we have presented an improved event- 
plane method for the flow analysis, which automatically 
corrects for autocorrelations and nonflow effects. As in 
the standard method, each event has its event plane ^P^, 
an estimate of the reaction plane, which is the same as for 
the standard method, except for technical details in the 
practical implementation. The trick which removes auto- 
correlations and nonflow effects is that there is in addition 
an event weight. Anisotropic flow w„ is then estimated as 
a weighted average of cos(n(0 — A straightforward 

application of this method would be to measure jet pro- 
duction with respect to the reaction plane at LHC. With 
the traditional event-plane method, such a measurement 
would require to subtract particles belonging to the jet 



from the event plane; in addition, strong nonflow corre- 
lations are expected within a jet, which would bias the 
analysis. 
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APPENDIX A: PRACTICAL IMPLEMENTATION 

Before we describe the implementation of the method, 
let us mention that there are in fact two Lee- Yang-zeroes 
methods, depending on how the generating function is 
deflned: the "sum generating function" makes explicit 
use of the flow vector [18| , while the "product generating 
function" [s^l is constructed using the azimuthal angles 
of individual particles, and cannot be expressed simply 
in terms of the flow vector. Cumulants also exist in both 
versions, the "sum" [l3| and the "product" [ll[. For 
Lee- Yang zeroes, both the sum and the product give es- 
sentially the same result for the lowest harmonic [15] : the 
difference between results from the two methods is signif- 
icantly less than the statistical error. On the other hand, 
the product generating function is significantly better 
than the sum generating function if one analyzes V4 or 
vi [26| using mixed harmonics. The method described 
below is strictly equivalent to the sum generating func- 
tion, although expressed in different terms. On the other 
hand, the product generating function cannot be recast 
in a form similar to the event-plane method, and will not 
be used here. 

The method a priori requires two passes through the 
data, which are described in Sec. lA II and Sec. IA2I 

1. First pass: locating the zeroes 

As with other flow analyses, one must first select events 
in some centrality class. The whole procedure described 
below must be carried out independently for each cen- 
trality class. 

The flow vector {Qx,Qy) is defined by Eq. (g]). In 
contrast to the standard event-plane method, no flatten- 
ing procedure is required to make the distribution of Q 
isotropic. Corrections for azimuthal anisotropies in the 
acceptance, which do not vary significantly in the event 
sample used, are handled using the procedure described 
below. We do not define the event plane as the az- 
imuthal angle of the flow vector, as in Eq. The pro- 
cedure below deflnes both the event weight and the event 
plane. 

The analysis uses the projection of the flow vector 
onto an arbitrary direction, see Eq. In practice, 

the first pass should be repeated for several equally- 
spaced values of n9 between and tt. This reduces 
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the statistical error, 5 values of 6 suffice, see Eq. IA5I 
For elliptic flow, for instance, 9 takes the values 9 = 
0,7r/10,27r/10,37r/10,47r/10. 

One first computes the modulus |G^(r)|, with de- 
fined by Eq. ([7]), as a function of r for positive r. One 
determines numerically the first minimum of this func- 
tion. This is the Lee- Yang zero. We denote its value by 
r^. It must be stored for each 9. 



2. Second pass: determining the event weight, wr, 
and the event plane, "ifR. 

In the second pass, one computes and stores, for each 
9, the following complex number: 



(Al) 



where joi — 2.40483. Except for statistical fiuctuations 
and asymmetries in the detector acceptance, should 
be purely imaginary. 

For each event, the event weight and the event plane 
are defined by 



Wrcosu'^/r = ( Re 



W^ijsinn^'fl = ( Re 




sinn6' ) , (A2) 



where Re denotes the real part, and angular brackets 
denote averages over the values of 9 defined in subsection 
lA 11 Our estimate of w„, denoted by i;„{LYZ}, is then 
defined by Eq. ©. 

We now discuss how the angle defined by Eq. (|A2|) 
compares with the event-plane from the standard analy- 
sis. First, we note that Eqs. (|A2p uniquely determine the 
angle n^'^ (modulo 2tt) only if the sign of Wu is known. 
The simplest convention is Wji > 0. In the simpli- 
fied implementation described in Sec. |TT1 however, where 
coincides with the standard event plane, Wr defined 
by Eq. ^ can be negative, because the Bessel function 
changes sign (see Fig. [1]). The convention Wu > then 
leads to a value of n'i/u which differs from the standard 
event plane by tt, since changing the sign of Wr amounts 
to shifting ti'^/r by tt in Eqs. (|A2[) . This is illustrated in 
Fig. [21 which shows the distribution of the relative angle 
between 4* r and the standard event plane in the simula- 
tion of V2 at LHC described in Sec. IIIII The distribution 
has two sharp peaks at and 7r/2. The sign ambiguity 
produces the peak at tt/2. The width of the peaks results 
from statistical fluctuations. The final result u„{LYZ}, 
given by Eq. ([3]), docs not depend on the sign chosen for 
Wh. 

If one wishes to have an event-plane as close as possible 
to the standard event plane, one may choose the following 
convention. Denoting by 5'^'^ the standard event plane. 



in 

o 
u 



10^ 



V„ = 0.0625 

— X = l-5 

— X = l- 



0.5 



1.5 
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, [rad] 



FIG. 3: Distribution of the relative angle between the event 
plane ^'h defined by Eq. (|A2|I . with Wr > 0, and the standard 
event plane, for the reconstruction shown in Fig. [J] 



one computes the following quantity: 

S = WB.cosn^RCOsn^%'^ + Wr sin n'^ r sin , 

(A3) 

where Wrcosti'^r and WRsinn'^R are defined by 
Eq. (|A2[) . The sign of Wr is then chosen as the sign 
of S, which ensures that n^R — n^E'^'^ lies between — 7r/2 
and 7r/2. 

The procedure described in this Appendix differs from 
the procedure described in Sec. [Til only in the case of 
non-uniform acceptance. This agreement can be seen 
in Fig. [1] which displays a comparison between the 
two. The solid line corresponds to the weight defined 
in Sec.[TTl(Eqs. ^ and (|12p). while the stars corresponds 
to the weight defined by Eq. (|A2p . as implemented in 
the Monte-Carlo simulation presented in Sec. Illll The 
agreement is very good. This agreement can also be seen 
directly on the equations. If the detector has perfect 
azimuthal symmetry, and in Eq. (|A2p arc inde- 
pendent of 9, up to statistical fluctuations. Neglecting 
these fluctuations, replacing with Eq. (O and inte- 
grating over 9, one easily recovers Eq. If there are 
azimuthal asymmetries in the detector acceptance, on 
the other hand, they are automatically taken care of by 
Eq. (|A2[) . The fact that one first projects the flow vector 
onto a fixed direction 9 is essential (for a related discus- 
sion, see [31II). 



3. Statistical errors 

The statistical error strongly depends on the resolution 
parameter 9] %, which is closely related to the reaction 
plane resolution in the event-plane analysis. It is given 
by 



X ■ 



(A4) 



In this equation, Vn is given by Eq. ([5]) , averaged over 9 
to minimize the statistical dispersion. The average val- 
ues (Qx), (Qj/), (Qx) ^iid {QD must be computed in the 
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first pass through the data. Please note that {Qx) and 
{Qy) vanish for a symmetric detector: they are accep- 
tance corrections. 

The price to pay for the ehmination of nonfiow effects 
is an increased statistical error. This increase is very 
modest if x is larger than 1: If x = 1-5, the error is only 
25% larger than with the standard event-plane method. 
If X = li it is larger by a factor 2. If x = 0.6, it is 20 
times larger. This prevents the application of Lee- Yang 
zeroes in practice for x smaller than 0.6. 

We now recall the formulas which determine the 
statistical error on w„{LYZ}: 



1 



4iV'Ji(joi)2p^^' 



kn 
P 



exp(^^cos 



exp 



2x^ 



■ cos 



fc7r\ 
P 
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(A5) 



where N' denotes the number of of objects one correlates 
to the event plane, whatever they are (jets, individual 
particles), and p is the number of equally-spaced values 
of 9 used in the analysis (see above). The larger p, the 
smaller the error. The recommended value is p = 5, be- 
cause larger values do not significantly reduce the error. 
This equation shows that the statistical error diverges 
exponentially when x is small. 
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